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We discuss data analysis techniques that can be used in the search for gravitational wave bursts 
from cosmic strings. When data from multiple interferometers are available, we describe consis- 
tency checks that can be used to greatly reduce the false alarm rates. We construct an expression 
for the rate of bursts for arbitrary cosmic string loop distributions and apply it to simple known 
solutions. The cosmology is solved exactly and includes the effects of a late-time acceleration. We 
find substantially lower burst rates than previous estimates suggest and explain the disagreement. 
Initial LIGO is unlikely to detect field theoretic cosmic strings with the usual loop sizes, though it 
may detect cosmic superstrings as well as cosmic strings and superstrings with non-standard loop 
sizes (which may be more realistic). In the absence of a detection, we show how to set upper limits 
based on the loudest event. Using Initial LIGO sensitivity curves, we show that these upper limits 
may result in interesting constraints on the parameter space of theories that lead to the production 
of cosmic strings. 
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I. INTRODUCTION 

Cosmic strings can form during phase transitions in 
the early universe Although cosmic microwave back- 
ground data has ruled out cosmic strings as the primary 
source of density perturbations, they are still potential 
candidates for the generation of a host of interesting as- 
trophysical phenomena including gravitational waves, ul- 
tra high energy cosmic rays, and gamma ray bursts. For 
a review see |3 . 

Following formation, a string network in an expand- 
ing universe evolves toward an attractor solution called 
the "scaling regime" (see "^l and references therein). In 
this regime, the statistical properties of the system, such 
as the correlation lengths of long strings, scale with the 
cosmic time, and the energy density of the string net- 
work becomes a small constant fraction of the radiation 
or matter density. 

The attractor solution arises due to reconnections, 
which for field theoretic strings essentially always oc- 
cur when two string segments meet. Reconnections lead 
to the production of cosmic string loops, which in turn 
shrink and decay by radiating gravitationally. This takes 
energy out of the string network, converting it to grav- 
itational waves. If the density of strings in the network 
becomes large, then strings will meet more often, thereby 
producing extra loops. The loops then decay gravitation- 
ally, removing the surplus energy from the network. If, 
on the other hand, the density of strings becomes too 
low, strings will not meet often enough to produce loops, 
and their density will start to grow. Thus, the network 
is driven toward a stable equilibrium. 

Recently, it was realised that cosmic strings could 



be produced in string-theory inspired inflation scenarios 
H, ^ IE IBi Q ' These strings have been dubbed cosmic 
superstrings. The main differences between field theo- 
retic strings and superstrings are that the latter recon- 
nect when they meet with probabilities p that can be less 
than 1, and that more than one kind of string can form. 
The values suggested for p are in the range 10"'^ — 1 j^]- 
The suppressed reconnection probabilities arise from two 
separate effects. The first is that fundamental strings in- 
teract probabilistically, and the second that strings mov- 
ing in more than 3 dimensions may more readily avoid 
reconnections: Even though two strings appear to meet 
in 3 dimensions, they miss each other in the extra dimen- 
sions 0. 

The effect of a small reconnection probability is to in- 
crease the time it takes the network to reach equilibrium, 
and to increase the density of strings p at equilibrium. 
The amount by which the density is increased is the sub- 
ject of debate, wit h p redictions ranging from p oc 
to p c>c p^"-^ |y, In this work we will assume the 

density scales like p oc p~^, and that only one kind of 
cosmic superstring is present 11]. 

Cosmic strings and superstrings can produce powerful 
bursts of gravitational waves. These bursts may be de- 
tectable by first generation interferometric earth-based 
gravitational wave detectors such as LIGO and VIRGO 
[ill IT^ IT^. Thus, the exciting possibility arises that a 
certain class of string theories may have observable con- 
sequences in the near future [lj|. It should be noted 
that if the hierarchy problem is solved by supersymme- 
try then the detectability of bursts produced by cosmic 
superstrings is strongly constrained due to dilaton emis- 
sion 0. 
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The bursts we are most likely to be able to detect 
are produced at cosmic string cusps. Cusps are regions 
of string that acquire phenomenal Lorentz boosts, with 
gamma factors 7^1. The combination of the large 
mass per unit length of cosmic strings and the Lorentz 
boost results in a powerful burst of gravitational waves. 
The formation of cusps on cosmic string loops and long 
strings is generic, and their gravitational waveforms are 
simple and robust 0|. 

The three LIGO interferometers (IFOs) are operating 
at design sensitivity. The VIRGO IFO is in the final 
stages of commissioning. Thus, there arises a need for a 
quantitative and observationally driven approach to the 
problem of bursts from cosmic strings. This paper at- 
tempts to provide such a framework. 

We consider the data analysis infrastructure that can 
be used in the search for gravitational wave bursts from 
cosmic string cusps, as well as cosmological burst rate 
predictions. 

The optimal method to use in the search for signals 
of known form is matched-filtering. In Sec. II we dis- 
cuss the statistical properties of the matched-filter out- 
put, template bank construction, an efficient algorithm 
to compute the matched-filter, and, when data from mul- 
tiple interferometers is available, consistency checks that 
can be used to greatly reduce the false alarm rate. In 
Sec. Ill we discuss the application of the loudest event 
method to set upper limits on the rate, the related issue 
of sensitivity, and background estimation. In Sec. IV we 
derive an expression for the rate of bursts for arbitrary 
cosmic string loop distributio ns, and app ly it to the dis- 
tribution proposed in Refs. 0, hj. In Sec. V we 
compare our estimates for the rate with results suggested 
by previous estimates (in 0,0,0). We find substan- 
tially substantially lower event rates. The discrepancy 
arises primarily from our estimate of a detectable burst 
amplitude, and the use of a cosmology that includes the 
late-time acceleration. We also compute the rate for sim- 
ple loop distributions, using the results of Sec. IV and 
show how the detectability of bursts is sensitive to the 
loop distribution. In Sec. VI we discuss how, in the ab- 
sence of a detection, it is possible to use the upper limits 
on the rate discussed in Sec. Ill to place interesting con- 
straints on the parameter space of theories that lead to 
the production of cosmic strings. We show an example 
of these constraints. We summarise our results and con- 
clude in Sec. VII. 



II. DATA ANALYSIS 



The low frequency cutoff of the gravitational wave signal, 
/; , is determined by the size of the feature that produces 
the cusp. Typically this scale is cosmological. In prac- 
tice, however, the low frequency cutoff of detectable ra- 
diation is determined by the low frequency behaviour of 
the instrument (for the LIGO instruments, for instance, 
by seismic noise) . The high frequency cutoff, on the other 
hand, depends on the angle 9 between the line of sight 
and the direction of the cusp. It is given by. 



fh 
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(2) 



where L is the size of the feature that produces the cusp, 
and in principle can be arbitrarily large. The amplitude 
parameter of the cusp waveform A is. 



2/3 



(3) 



where G is Newton's constant, /i is the mass per unit 
length of strings, and r is the distance between the cusp 
and the point of observation. We have taken the speed 
of light c = 1 . 

The time domain waveform can be computed by taking 
the inverse Fourier transform of Eq. . The integral has 
a solution in terms of incomplete Gamma functions. For 
a cusp with peak arrival time at t = it can be written 



h{t) =2T:A\t\^/^ { iT{~l/3,2TTifit)+c.c. 

- [iri-l/3,27rifht)+c.c.]}. (4) 

The duration of the burst is on the order of the inverse 
of the the low frequency cutoff // , and the spike at t = 
is rounded on a timescale ~ 

Eq. I0J is useful to make a rough estimate the am- 
plitude of changes in the strain h{t) that a cusp would 
produce. If we consider a time near t = 0, we can expand 
the incomplete Gamma functions according to 17j, 



T{a, x) — T{a) 



(-1)"2;"+" 



n\{a + n) 



(5) 



Keeping the first term of the sum we see that the ampli- 
tude in strain of the cusp is. 



« 6Af-'/\ h » fi. 



(6) 



A. The signal 



B. Definitions and conventions 



In the frequency domain, the waveforms for bursts of 
gravitational radiation from cosmic string cusps are given 

by El, 



The optimal method to use in the search for signals of 
known form is matched-filtering. For our templates we 
take 



h{f) = A\f\-^'-'Q{h^ J)Q{f ~ fi). 



(1) 



t(/) = /~'/'e(A-/)e(/-/o, 



(7) 
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so that a signal would have h(f) = At{f). We define the 
usual detector- noise- weighted inner product [2^ . 



(8) 



where Sh{f) is the single-sided spectral density defined 
by (n{f)n*if')) = ^Sif - nShif), where n(/) is the 
Fourier transform of the detector noise. 

The templates of Eq. 10 can be normalised by the 
inner product of a template with itself, 



so that, 



a' = m 



i=t/a, and 



(9) 



(10) 



If the output of the instrument is s(i), the signal to 
noise ratio (SNR) is defined as, 



p={s\i). 



(11) 



In general, the instrument output is a burst h{t) plus 
some noise n(i), s{t) = h{t) -f n{t). When the signal is 
absent, h{t) = 0, the SNR is Gaussian distributed with 
zero mean and unit variance, 

(p) = {{n\t)) = 0, {p") = {{n\tf) = 1. (12) 

When a signal is present, the average SNR is 

{p) = {{h\i)) + {{n\i)) = {Aai\t) = An, (13) 
and the fluctuations have variance 

{p')-{p?^{{h + n\if)-A^a^^l. (14) 

Thus, when a signal of amplitude A is present, the sig- 
nal to noise ratio measured, p , is a Gaussian random 
variable, with mean Aa and unit variance. 



p = Acj±\ 



(15) 



The amplitude that we assign the event A depends on 
the template normalisation cr. 



7 j5 . 1 
A = - = A± -, 



(16) 



This means that in the presence of Gaussian noise the 
relative difference between measured and actual signal 
amplitudes will be proportional to the inverse of the SNR, 



AA _ 1 



(17) 



If an SNR threshold pth is chosen for the search, on av- 
erage only events with amplitude 



(18) 



will be detected. 



C. Template bank 

Cosmic string cusp waveforms have two free parame- 
ters, the amplitude and the high frequency cutoff. Since 
the amplitude is an overall scale, the template need not 
match the signal amplitude. However, the high frequency 
cutoff does affect the signal morphology, so a bank of 
templates is needed to match possible signals. 

The template bank, {ii}, where i = l,2,...iV, is one- 
dimensional, and fully specified by a collection of high fre- 
quency cutoffs {fi}. The number of templates N should 
be large enough to cover the parameter space finely. We 
choose the ordering of the templates so that fi > fi+i- 

Although in principle the high frequency cutoff can 
be infinite, in practice the largest distinguishable high 
frequency cutoff is the Nyquist frequency of the time- 
series containing the data, /Nyq = 1/(2 At), where At is 
the sampling time of that time-series. The normalisation 
factor for the first template is given by 



{ti\ti 



/Nyq /-8/3 



(19) 



This template is the one with the largest cr, and thus the 
largest possible SNR at fixed amplitude. 

To determine the high frequency cutoffs of the remain- 
ing templates, we proceed as follows. The fitting factor 
between two adjacent templates ti and ti+-\ , with high 
frequency cutoffs /; and fi+i, is defined as [23] 



F 



{u\u 



\/ {U\ti){ti+l\ti+l) 



1 



(20) 



where e is the maximum mismatch we choose to allow. 
The maximum mismatch is twice the maximum fractional 
SNR loss due to mismatch between a template in the 
bank and a signal, and should be small. Since ti+i has a 
lower frequency cutoff than ti , 



{ti\ti+l) — {ti+l\ti+l)- 

Thus, the maximum mismatch is. 



1 



{u\u 



1 



(21) 



(22) 



and once fixed, determines the high frequency cutoff /i+i, 
given fi. Thus the template bank can be constructed 
iteratively. 



D. Single interferometer data analysis 

A convenient way to apply the matched-filter is to 
use the inverse spectrum truncation procedure p3 |. The 
usual implementation of the procedure involves the cre- 
ation of an FIR (Finite Impulse Response) digital filter 
for the inverse of the single-sided power spectral density. 
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Sh{f), of some segment of data. This filter, along with 
the template can then be efficiently applied to the data 
via FFT (Fast Fourier Transform) convolution. Here we 
describe a variation of this method in which the template 
as well as Sh{f) are incorporated into the digital filter. 

Digital filters may be constructed for each template j 
in the bank as follows. First we define a quantity which 
corresponds to the square root of our normalised tem- 
plate divided by the single-sided power spectral density, 

This quantity is then inverse Fourier transformed to cre- 
ate, 

Tj,i/2it) - J dfe'-'f%,,/,if). (24) 

At this point the duration of the filter is identical to 
the length of data used to estimate the power spectrum 
Sh{f)- We then truncate the filter, setting to zero values 
of the filter sufficiently far from the peak. This is done 
to limit the amount of data corrupted by filter initialisa- 
tion. After truncation, however, the filter still needs to 
be sufficiently long to adequately suppress line noise in 
the data (such as power line harmonics). The truncated 
filter is then forward Fourier transformed, and squared, 
creating the frequency domain FIR filter Tj (/) . The SNR 
is the time-series given by, 

/>OC 

P,(0=45R/ df s{f)T,{f)e^^^f\ (25) 
Jo 

where s(/) is a Fourier transform of the detector strain 
data. 

To generate a discrete set of events or triggers we can 
take the SNR time series for each template Pj{t) and 
search for clusters of values above some threshold pth- 
We identify these clusters as triggers and determine the 
following properties: 1) The SNR of the trigger which is 
the maximum SNR of the cluster; 2) the amplitude of 
the trigger, which is the SNR of the trigger divided by 
the template normalisation <Tj] 3) the peak time of the 
trigger, which is the location in time of the maximum 
SNR; 4) the start time of the trigger, which is the time 
of the first SNR value above threshold; 5) the duration of 
the trigger, which includes all values above the threshold; 
and 6) the high frequency cutoff of the template. The 
final list of triggers for some segment of data can then 
be created by choosing the trigger with the largest SNR 
when several templates contain triggers at the same time 
(within the durations of the triggers). 

The rate of events we expect can be estimated for the 
case of white Gaussian noise. In the absence of a sig- 
nal the SNR is Gaussian distributed with unit variance. 
Therefore, at an SNR threshold of pth, the rate of events 
with p > pth we expect with a single template is 

R{> pth) ~ At-iErfc(pth/V2) (26) 



where Ai is the sampling time of the SNR time series, and 
Erfc is the complementary error function. For example, 
if the time series is sampled at 4096 Hz and an SNR 
threshold of 4 is used we expect a false alarm rate of R{> 
4) K, 0.26 Hz. The rate of false alarms rises exponentially 
with lower thresholds. 



E. Multiple interferometer consistency checks 

When data from multiple interferometers is available, 
a coincidence scheme can be used to greatly reduce the 
false alarm rate of events. For example, in the case of 
the three LIGO interferometers, using coincidence would 
result in a coincident false alarm rate of 

i? = i?Hii?H2i?Li(2ArHH)(2ArHL), (27) 

where i?Hi and i?H2 are the rate of false triggers in the 
4 km and 2 km IFOs in Hanford, WA, and _Rli is the rate 
of false triggers in the 4 km IFO in Livingston, LA. The 
two time windows Athh, and Athl are the maximum al- 
lowed time differences between peak times for coincident 
events in the HI and H2 IFOs, and HI-LI (or H2-LI), 
respectively. Eq. H27|l follows trivially if the events are 
Poisson distributed. The factors of 2 multiplying the 
time windows arise from the fact that an event in the 
trigger set of one instrument will survive the coincidence 
if there is an event in the other instrument within a time 
Ar on either side. The time windows should allow for 
light-travel time between sites, shifting of measured peak 
location due to noise, timing errors and calibration errors. 
For example, if we take Athl — 12 ms and Athh = 2 ms, 
which allows for the light travel time of 10 ms between the 
Hanford and Livingston sites as well as 2 ms for other un- 
certainties, and the single IFO false alarm rate estimated 
in the previous subsection (0.26 Hz), we obtain a triple 
coincident false alarm rate, 

i? « 1.7 X lO^'^Hz, (28) 

which is about 50 false events per year. 

Additionally, when two instruments are co-aligned, 
such as the HI and H2 LIGO IFOs, we can demand some 
degree of consistency in the amplitudes recovered from 
the matched-filter (see, for example, the distance cut used 
in ^3)- If instruments are not co-aligned then the 
amplitudes recovered by the matched-filter could be dif- 
ferent. This is due to the fact that the antenna pattern 
of the instrument depends on the source location in the 
sky. If the event is due to a fluctuation of the noise or an 
instrumental glitch wc do not generally expect the am- 
plitudes in two co-aligned IFOs to agree. The effect of 
Gaussian noise to I cr on the estimate of the amplitude 
can be read from Eq. (|I7|I . This equation can be used to 
construct an amplitude veto. For example, to veto events 
in HI we can demand 

2 < ± + « h (29) 

^Hl VPHl / 
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where S is the number of standard deviations we allow 
and K is an additional fractional difference in the ampli- 
tude that accounts for other sources of uncertainty (such 
as the calibration). 

These tests can be used to construct the final trigger 
set, a list of survivors for each of our instruments which 
have passed all the consistency checks. 



III. UPPER LIMITS, SENSITIVITY AND 
BACKGROUND 

In this section we describe the application of the loud- 
est event method [2^ to our problem, address the related 
issue of sensitivity, and discuss background estimation. 

We define the loudest event to be the survivor with the 
largest amplitude. Suppose, in one of the instruments, 
we find the loudest event to have an amplitude A*. For 
an observation time T, on average, the number of events 
we expect to show up with an amplitude greater than A* 



where the effective rate is defined as 
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eiA*,A)'-^dA. 



(30) 



(31) 



Here, e{A* , A) is the search efficiency, namely the fraction 
of events with an optimally oriented amplitude A found 
at the end of the pipeline with an amplitude greater than 
A*. Aside from sky location effects, the measured am- 
plitude of events is changed as a result of the noise (see 
Eq. HI?!) ). The efficiency can be measured by adding sim- 
ulated signals to detector noise, and then searching for 
them. The efficiency accounts for the properties of the 
population, such as the distribution of sources in the sky 
and the different high frequency cutoffs (as we will see, 
dR/df cx f^^^^ for cusps), as well as possible inefficien- 
cies of the pipeline. The quantity 



dRjA) 
dA 



dA 



(32) 



is the cosmological rate of events with optimally oriented 
amplitudes between A and A+dA and will be the subject 
of the next Section. 

If the population produces Poisson distributed events, 
the probability that no events show up with a measured 
amplitude greater than A* , in an observation time T, is 



P 



-T7 



(33) 



The probability, 77, that at least one event with amplitude 
greater than A* shows up is 77 = 1 — P, so that if = 0.9, 
90% of the time we would have expected to see an event 
with amplitude greater than A*. 
From Eq. we see that. 



ln(l - 77) -T7, 



(34) 



Setting 1] = 0.9, say, we have 



790% 



2.303 

T 



(35) 



We say that the value of 7 < 799% with 90% confidence, 
in the sense that if 7 = 790%, 90% of the time we would 
have observed at least one event with amplitude greater 
than A*. If we take a nominal observation time of one 
year, T = 3.2 x 10^ s, our upper limit statement becomes. 



7 < 790% ~ 2.303/year « 7.3 x 10~®s~ 



(36) 



This bound on the effective rate can be used to con- 
strain the parameters of cosmic string models that enter 
through the cosmological rate dR/dA in Eg. 1311 

The efficiency e{A*,A) is the fraction of events with 
optimally oriented amplitude A which show up in our 
final trigger set with an amplitude greater than A*. If 
A:$> A* then we expect e(^*, A) « 1, and when A < A*, 
e{A*,A) Ri 0. If ^50% is the amplitude at which half 
of the signals are found with an amplitude greater than 
A*, we can approximate the sigmoid we expect for the 
efficiency curve by. 



eiA\A)^eiA-A,o%), 



(37) 



meaning we assume all events in the data with an am- 
plitude A > AcjQ% survive the pipeline and are found 
with an amplitude greater than A*. This means we can 
approximate the rate integral 7 of Eq. H31|l as. 
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dR{A) 
dA 



dA^ R 



(38) 



so that if R:^A^o% ^ 790% for a particular model of cos- 
mic strings, the model can be excluded (in the frequentist 
sense) at the 90% level. We expect ^50% to be propor- 
tional to A* . It is at least a factor of ^/5 larger, due to the 
averaging of source locations on the sky [23 , and poten- 
tially even larger than this due inefficiencies of the data 
analysis pipeline. If there are no significant inefficiencies 
then we expect ^50% ~ ^/EA* . 

At design sensitivity the form of the Initial LIGO noise 
curve may be approximated by [28| . 



Shif) = 



28 



1.09 X 10" 



1.44 X 10 



-45 



^30Hz^ 



+ 1.28 X lO"''^ 1 



/ 



90Hz, 



s. (39) 



The first term is the so-called seismic wall, the second 
term comes from thermal noise in the suspension of the 
optics, and the third term from photon shot noise. In the 
absence of a signal, this curve can be used to estimate 
the amplitude of typical noise induced events. 
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As pointed out in i a reasonable operating point for 
the pipeline involves an SNR threshold of pth = 4. The 
value of the amplitude is related to p, the SNR, and u, 
the template normalisation, 

'^"i *^^^/. ''km <^°' 

by Eq. ((TH|l . If the loudest surviving event has an SNR 
close to the threshold, ~ 4, a low frequency cutoff at 
/; — 40 Hz (the seismic wall is such that it's not worth 
running the matched- filter below 40 Hz), and a typical 
high frequency cutoff fh = 150 Hz then the amplitude 
of the loudest event wiU be A» w 4 x lO^^^s^^/^. So we 
expect ^50% « 9 X IQ-^^s"^/^ for the Initial LIGO noise 
curve. 

The question of the sensitivity of the search is subtly 
different. In this case we are interested in how many 
burst events we can detect at all, not just those with 
amplitudes greater than A* . The rate of events we expect 
to be able to detect is thus, 

^s^l^e{0,A)^dA. (41) 

The efhciency includes events of all amplitudes and if the 
cosmological rate dR/dA favours low amplitude events, 
this rate could be larger than that given by Eq. (|30|l . 
Furthermore this rate should be compared with a rate 
of 1/T, rather than, say, the 90% rate of Eq. The 
net effect is that the parameter space of detectability 
can be substantially larger than the parameter space over 
which we can set upper limits. It should be noted that 
the in the fraction of detected events e(0. A) in Eq. HTl 
means that we only demand that events with optimally 
oriented amplitudes A show up in the final trigger (not 
that their amplitudes be larger than A*). There is in 
fact a minimum detectable amplitude of events, given 
by A,„in = pth/o-max, where CTmax is the largest template 
normalisation used in the search. 

In the remainder of the paper we will use an amplitude 
estimate of ^50% = IQ-'^^s-^'^ for Initial LIGO upper 
limit and sensitivity estimates. Advanced LIGO is ex- 
pected to be somewhat more than an order of magnitude 
more sensitive then Initial LIGO, so for the Advanced 
LIGO rate estimates we will use ^50% = 10~^^s~^/'^. It 
should be noted that the loudest event resulting from a 
search on IFO data could be larger than this, for instance 
if it was due to a real event or a surviving instrumental 
glitch. 

The rate of events that we expect in the absence of 
a signal, called the background, can be estimated either 
analytically (via Eq. 1)27(1 . in which the single IFO rates 
are measured), or by performing time-slides (see, for ex- 
ample )■ In the latter procedure single IFO triggers 
are time-shifted relative to each other and the rate of ac- 
cidental coincidences is measured. Care must be taken 
to ensure the time-shifts are longer than the duration of 
real events. If the rate of events is unknown, then one 



can look for statistically significant excesses in the fore- 
ground data set. In this case a statistically significant 
excess could point to a detection. However, as we shall 
see, there is enough freedom in models of string evolu- 
tion to lead to statistically insignificant excesses in any 
background. Therefore it is prudent to carefully exam- 
ine the largest amplitude surviving triggers regardless of 
consistency with the background. 

IV. THE RATE OF BURSTS 

In this section we will derive the expression for the cos- 
mological rate as a function of the amplitude, which is 
needed to evaluate the upper limit and sensitivity inte- 
grals (Eqs. H30|l and H41|) '). To make this account self- 
contained we will re-derive a number of the results pre- 
sented in 0, [T3 . IT^ . The main differences from the 
previous approach are that we derive an expression for a 
general cosmic string loop distribution that can be used 
when such a distribution becomes known, we absorb some 
uncertainties of the model, as well as factors of order 
1, into a set of ignorance constants, and we consider a 
generic cosmology with a view to including the effects of 
a late time acceleration. 

Numerical simulations show that at any cosmic time 
i, the network consists of a few long strings that stretch 
across the horizon and a large collection of small loops. 
In early simulations 18] the size of loops, as well as the 
substructure on the long strings, was at the size of the 
simulation resolution. At the time of the early simu- 
lations, the consensus reached was that the size of loops 
and sub-structure on the long strings was at the size given 
by grav itational back-reaction. Very recent simulations 
[19I I20I I21I I suggest that the size of loops is given by 
the large scale dynamics of the network, and is thus un- 
related to the gravitational back-reaction scale. Further 
work is necessary to establish which of the two possibili- 
ties is correct. 

In either case, the size of loops can be characterised by 
a probability distribution, 

n{l,t)dl, (42) 

the number density of loops with sizes between I and l+dl 
at a cosmic time t. This distribution is unknown, though 
analytic solutions can be found in simple cases (see [3| 
§9.3.3). 

The period of oscillation of a loop of length I isT ~ 1/2. 
If, on average, loops have c cusps per oscillation, the 
number of cusps produced per unit space-time volume 
by loops with lengths in the interval dl at time t is |35| 

i^{l,t)dl^^n{l,t)dl. (43) 
We can write the cosmic time as a function of the redshift, 

t = H^'M^), (44) 
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where Hq is current value of the Hubble parameter and 
ftiz) is a function of the redshift z which depends on the 
cosmology. In 0, 0, ^| , an approximate interpolating 
function for (ftiz) was used for a universe which contains 
matter and radiation. For the moment we will leave it 
unspecified. Later we will compute it numerically for 
a more realistic cosmology which includes the late-time 
acceleration. We can therefore write the number of cusps 
per unit space-time volume produced by loops with sizes 
between / and I -\- dl at a redshift z, 



2c 

i>{l,z)dl — —n{l, z)dl. 



(45) 



Now we would like to write the length of loops in terms 
of the amplitude of the event an optimally oriented cusp 
from such a loop would produce at an earth-based de- 
tector. The amplitude we expect from a cusp produced 
at a redshift z, from a loop of length I can be read off 
Eq. (4.10) of P|. It is 



A 
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(H-z)i/3^^(z)- 



(46) 



Here gi is an ignorance constant that absorbs the uncer- 
tainty on exactly how much of the length I is involved in 
the production of the cusp, as well as factors of 0(1) that 
have been dropped from the calculation (see the deriva- 
tion leading up to Eq. (3.12) in 0). If loops are not 
too wiggly, we expect this ignorance constant to be of 
L). We have expressed r{z), the amplitude distance 
as an unspecified cosmology dependent function of 
the redshift ipr{z), 



The amplitude distance is a factor of (1 - 
the luminosity distance. Since, 



-z) smaller than 



l{A,z) 



3/2 



and — 7 ~ — -I, 
dA 2A ' 



(48) 

we can write the number of cusps per unit space-time 
volume which produce events with amplitudes between 
A and A + dA from a redshift z as, 



3c 

i^{A, z)dA = —n{l{A, z), z)dA. 



(49) 



As discussed in IT2| , we can observe only a fraction of 
all the cusps that occur in a Hubble volume. Suppose 
we look, via matched-filtering, for signals of the form of 
Eq. (Q). The highest frequency we observe fh is related 
to the angle 9 that the line of sight makes to the direction 
of the cusp. If /* is the lowest high frequency cutoff that 
our instrument can detect with confidence, the maximum 
angle the line of sight and the direction of a detectable 
cusp can subtend is, 



(52/*0 



-1/3 



(50) 



The ignorance constant 52 absorbs a factor of about 2.31 
[l3l |. other factors of 0(1), as well as the fraction of the 
loop length / that actually contributes to the cusp (see the 
derivation leading up to Eq. (3.22) of H^). Just like for 
our first ignorance constant gi , if loops are not too wiggly, 
we expect 32 to be of 0{1). For bursts originating at 
cosmological distances, red-shifting of the high frequency 
cutoff must be taken into account. So we write, 



,(z,^,^) = [52(l + ^)/*r'/^ 



(51) 



The waveform in Eq. Jy), was derived in the limit 9 <C 
1, so that cusps with 61 > 1 will not have the form of 
Eq. (indeed, they are not bursts at all). Therefore 
cusps with 6' > 1 should not contribute to the rate. Thus, 
the fraction of cusps that are observable, A, is 



eil-9rn{zj,,l)). (52) 



The first term on the right hand side is the beaming 
fraction for the angle 9m- At fixed /, the beaming fraction 
is the fraction of cusp events with high frequency cutoffs 
greater than The 8 function cuts off cusp events 
that have 9m > 1; it ensures we only count cusps whose 
waveform is indeed given by Eq. 

The rate of cusp events we expect to see from a volume 
dV{z) (the proper volume in the redshift interval dz), in 
an amplitude interval dA is |36j 



dR 



dV{z)dA 



{l + z)-'u{A,z)A{zJ,,A), 



(53) 



where A{z,f^,A) = A{z, f^,l{A, z)). The factor of 
(47) (1 + z)^^ comes from the relation between the observed 



burst rate and the cosmic time: Bursts coming from large 
redshifts are spaced further apart in time. We write the 
proper volume clement as. 



dV{z) 



Hq ^ipv{z)dz, 



(54) 



where (pv{z) is a cosmology dependent function of the 
redshift. Thus, for arbitrary loop distributions, we can 
write the rate of events in the amplitude interval dA, 
needed to evaluate the upper limit integral Eq. (PHI), and 
sensitivity integral Eq. (|41|) as. 



dA-"' 



, dzipv{z){l'V z) V(A,z)A(z,/,,^). 
Jo 

(55) 

Damour and Vilenkin 

[Ullllll took the loop distri- 
bution to be 



n(/,t) = {pTGii)-^t-^5{l - at), 



(56) 



where p is the reconnection probability, which can be 
smaller than 1 for cosmic superstrings. The constant F 
is the ratio of the power radiated into gravitational waves 
by loops to G/i^, and thus related to the lifetime of loops. 
It is measured in numerical simulations to be F 50. 
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In this loop distribution all loops present at a cosmic 
time t, are of the same size at. The distribution is con- 
sistent with the assumption usually made in the litera- 
ture that the size of loops is given by gravitational back- 
reaction, and that a ^ TGfi. Recently it was realised 
that damping of perturbations propagating on strings 
due to gravitational wave emission is not as efficient as 
previously thought jsj- As a consequence, the size of 
the small-scale structure is sensitive to the spectrum of 
perturbations present on the strings, and can be much 
smaller than the canonical value FG/it js^ • 

If the value of a is given by gravitational back-reaction 
we can write it as 

a^e {rGfiy . (57) 

We use two parameters e, introduced in as well as n 
that can be used to vary the size of loops. The parameter 
n arises naturally from gravitational back-reaction mod- 
els and is related to the power spectrum of perturbations 
on long strings [s^l- For example, if the spectrum is in- 
versely proportional to the mode number, then n = 3/2. 
This is the spectrum of perturbations we expect if, say, 
the shape of the string is dominated by the largest kink 

m 

As we have mentioned, recent simulations suggest that 
loops are produced at sizes unrelated to the gravitational 
back-reaction scale ^3113,1^. If this is the case, then 
the loops produced survive for a long time, and the form 
of the distribution can be computed analytically in the 
matter and radiation eras using some simple assumptions 
(see §9.3.3). 

In general the rate integral needs to be computed, pre- 
sumably numerically, through Eq. 155() . In the case of 
the more simple loop distribution of Eq. H56|) , it is conve- 
nient to proceed slightly differently. In this case all loops 
at a given redshift are of the same size, so we can directly 
associate amplitudes with redshifts. 

We can write the rate of cusp events we expect to see 
from the redshift interval dz, from loops of length in the 
interval dl as 



say Eq. H38(l . over redshifts rather than over amplitudes. 



dR 
dzdl 



H^\v{z){l + z)-^v{l, z)A{z, (58) 



Using Eq. (|45|1 . Eq. H5()|l with t expressed in terms of the 
redshift, and integrating over / we find. 



dR 

dz 



= Ha 



X 9(1 - 6l„j(z,/,,ai7o"Vt(^)) 



^r^^/^(z)^v(z)(l + z)-5/3 
(59) 



At a redshift of z, all cusps produce bursts of the 
same amplitude, given by replacing I with aHQ^ipt[z) 
in Eq. (|45|l . Therefore the solution of. 



2/3/ N 
"Pt (2) 



AH, 



'1/3 



(l + z)i/3^,(z) giG^a2/3 



(60) 



for z gives the redshift from which a burst of amplitude A 
originates. This means we can perform the rate integral. 



(61) 



where 259% is the solution of Eq. H60I) for A — A^q% . 

To compute the cosmological functions, ^pt{z), ^r{z), 
and ipv{z), we use the set of cosmological parameters 
in 34J, which provide a good fit to recent cosmological 
data. The precise values of the parameters are not crit- 
ical, but we include them here for clarity. They are, 
Hn ^ 73 km s^^ Mpc = 2.4 x 10'^^ s-\ n,^ = 0.25, 
ilr = 4.6x 10~^, and n\ = 1 — The derivation of 

the three cosmological functions is included in Appendix 
A for completeness and result in the cosmological func- 
tions Eqs. (|A4|) . HA6p . and HA8|I . In the remainder of the 
paper we will refer to this cosmological model as the A 
universe. 



V. RESULTS I: DETECTABILITY 

A. Comparison with previous estimates 

Damour and Vilenkin considered a uni- 
verse with matter and radiation. They introduced a set 
of approximate interpolating functions for ftiz), ipr{z), 
and ipv{z). They were 



iptiz) = toHoil + zy^'/^il + z/z,g)-'/^, 



1 + 



(62) 



(63) 



and 



Vv{z) = itoH„)HO^z^{l + z)-^^/^l + z/z,,)-^/\ (64) 



Damour and Vilenkin used z^q ~ lO"^'^ as the redshift of 
radiation matter equality and set the age of the universe 
tg = IQi''-^ s. Notice that if we substitute Eqs. (|^ and 
(|64|l into Eq. H59|l and set a = TGfi and g2 ~ p — 1 we 
recover Eq. (5.17) of |13| (aside from a factor of z related 
to their use of a logarithmic derivative). 

We will use Eqs. (ESI and to compare the 

previous results [llj, JJj with those of a cosmology 
solved exactly including the effects of a late time accel- 
eration, i.e. the cosmological functions Eqs. I)A4|I . ljA6|l . 
and (jA8|) . Damour and Vilenkin considered the ampli- 
tude (in strain) of events at a fixed rate of 1 per year and 
compared that to an SNR 1, optimally oriented event 
(see, for example, the dashed horizontal lines of Fig. 1 
in UM)- Instead, we have started by estimating a de- 
tectable signal amplitude, and computed the event rate 
at and above that amplitude. 

Bursts from cosmic string cusps are Poisson dis- 
tributed. The quantity we have computed, 7, is the rate 
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FIG. 1: Plot of the rate of gravitational wave bursts, 7 
(top panel), and the probability of having at least least 
one event in our data set with amplitude larger than ^50% 
in a year of observation (bottom panel), as a function of Gfi. 
For all curves we have set a = FG/i, F = 50, /, = 75 Hz, 
c — p = 1, and the ignorance constants gi =32 ~ 1- The 
dash-dot and dashed curves show 7 and rj computed with 
the Damour-Vilenkin cosmological functions Eqs. 1621 . I63II 
and with Aso% = 10"^^ s-^''^ and ^50% = 10"^° s"^/^ 
respectively. The thick and thin solid curves show 7 and 77 
computed in a universe with a cosmological constant with 
amplitudes ^50% = 10"^° s"^/^, and ^50% = 10"^^ s"^/^ 
respectively. 

of events in our data set with amplitudes greater than 
A^Q%. In an observation time T the probability of not 
having such an event is exp(— 7T). Hence, the odds of 
having at least least one event in our data set with am- 
plitude larger than our minimum detectable amplitude 
is, 

V = l~ e-^^. (65) 

Figure ^ shows the rate of burst events, 7, as well 
as the probability 77 of having at least least one event 
in our data set with amplitude larger than ^50% for a 
year of observation, as a function of Gfi for two different 
models. For all curves we have set a = FG/i, F = 50, 
/* — 75 Hz, c = p = 1, and the ignorance constants 
<7i = 52 = 1- We will refer to string models with these pa- 
rameters as "classic" , which is appropriate for field theo- 
retic strings with loops of size I = TG^t. The dashed-dot 



and dashed curves of Fig. ^ show 7 and rj computed us- 
ing the Damour-Vilenkin cosmological functions namely, 
Eqs. (EIJ), (ESI and For the dashed-dot curves we 

have used an amplitude estimate of ^50% — 10~^^ s~^/'^. 
This amplitude estimate can be obtained using the Ini- 
tial LIGO sensitivity curve, setting the SNR threshold to 
1, and assuming all cusp events are optimally oriented 
(as used for the dashed horizontal lines of Fig. 1 in T3.]). 
This is also our estimate for the amplitude in the case of 
Advanced LIGO. The dashed curves show 7 and r] com- 
puted with the amplitude A^q% = 10~^° s"^/"^, which we 
feel is more appropriate for Initial LIGO. The thick and 
thin solid curves show 7 and 77 computed by evaluating 
the cosmological functions (Eqs. (|A4p . IjAGp and IjASp l 
numerically for the A universe (see Appendix A). The 
thick solid curves correspond to our amplitude estimate 
for Initial LIGO, and the thin solid curves to our estimate 
for Advanced LIGO. 

The functional dependence of the rate of gravitational 
wave bursts on G/i is discussed in detail in Appendix B. 
Here we summarise those findings. From left to right, 
the first steep rise in the rates as a function of Gfi of 
the dashed and dashed-dot curves of Fig. ^ comes from 
events produced at small redshifts {z ^ 1). The peak and 
subsequent decrease in the rate starting around G/i ^ 
10~^ comes from events produced at larger redshifts but 
still in the matter era (1 <C z <C Zeq). The final rise comes 
from events produced in the radiation era {z ^ Zeq). 

For "classic" cosmic strings {p = e = n = 1 ) , the mat- 
ter era maximum in our estimate for the rate of events 
at Initial LIGO sensitivity is about 7 x 10~^ events per 
year, which is substantially lower than the rate ~ 1 per 
year su gges ted by the results of Damour and Vilenkin 
jlll IT^ The bulk of the difference arises from our 
estimate of a detectable amplitude. This is illustrated by 
the dashed-dot and dashed curves of Fig. ^ which use the 
same cosmological functions, and two estimates for the 
amplitude, ^50% = lO'^i s'^/^ and = 10~^° s'l/^ 

respectively. Our amplitude estimate results in a de- 
crease in the burst rate by about a factor of 100 at the 
matter era peak. A more detailed discussion of the effect 
of the amplitude on the rate can be found in Appendix B. 
The remaining discrepancy arises from differences in the 
cosmology, as well as factors of 0(1) that were dropped 
in the previous estimates, which account for a further 
decrease by factor of about 10. This is illustrated by 
the difference between the dashed-dot and thin solid 
curves of Fig. ^ which use the same amplitude estimate 
^50% = 10"^"'^ s~^/^, and the Damour-Vilenkin cosmo- 
logical interpolating functions (Eqs. Ht)2|) . Ht)3|) and (|64|) 'l 
and the A universe functions (Eqs. HA4|) . HA6|I and EU) 
respectively. When z << 1, the effects of a cosmological 
constant are un-important and differences arise from fac- 
tors of 0(1) that were dropped in the previous estimates. 
For z > 1, the differences arise from a combination of 
the effects of a cosmological constant as well as factors 
of 0(1). The net effect is that the chances of seeing an 
event from "classic" strings using Initial LIGO data have 
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FIG. 2: Comparison of 7 (top panel) and 77 with a year of 
observation (bottom panel), as a function of G/i for several 
string models. For all three curves an Initial LIGO ampli- 
tude estimate of A^q% = 10~^° s~^'^^ has been used, and the 
cosmological functions have been evaluated in the A universe. 
The model parameters are identical to those of Fig. except 
where indicated. The solid curves show 7 and rj computed 
with a reconnection probability of p = 10^^. The dashed 
curves show 7 and 77 computed with a size of loops given by 
Eq. (1571 with £ = 1, and n = 3/2. The dashed-dot curves 
show the combined effect of a low reconnection probability, 
p — 10~^, as well as a size of loops given by Eq. 1571 with 
£ = 1, and n — 3/2. 



dropped from order unity to about lO^'^ at the matter 
era peak. This is illustrated by the difference between the 
dashed-dot and thick solid curves of Fig. ^ The dashed- 
dot curves were computed using the amplitude estimate 
^50% = 10~^^ s~^/'^ and the Damour-Vilenkin inter- 
polating cosmological functions, whereas the thick solid 
curves use an amplitude estimate of ^50% — 10^^" s^^/^ 
in the A universe. 

Cosmic superstrings, however, may still be detectable 
by Initial LIGO. Furthermore, if the size of the small- 
scale structure is given by gravitational back-reaction, 
reasonable estimates for what the size of loops might 
be also lead to an enhanced rate of bursts. Figure [21 
illustrates this point. All curves use the Initial LIGO 
amplitude estimate of ^50% — 10^^° s^^/'^, and the cos- 
mological functions Eqs. (jA4|) . (jA6|) and ljA8|l computed 
in the A universe. The solid curves show 7 and rj com- 



FIG. 3: Same as Fig.|5|but with the Advanced LIGO ampli- 
tude estimate of ^59% = 10^^""^ s""'^''^. 



puted with a reconnection probability of p = 10~^. The 
dashed curves show 7 and i] computed for loops with a 
size given by Eq. H57|) with e = 1, and n — 3/2. This 
is the value of n we expect when the spectrum of per- 
turbations on long strings is inversely proportional to 
the mode number, and the result we expect if the spec- 
trum of perturbations on long strings is dominated by the 
largest kink [s^- The dashed-dot curves show the com- 
bined effect of a low reconnection probability, p = 10^"^, 
as well as a size of loops given by Eq. (|57|) with e = 1, 
and n — 3/2. The remaining parameters for all three 
curves are identical to those of Fig. ^ Advanced LIGO 
has a considerably larger chance of making a detection of 
cosmic superstrings or field-theoretic strings if loops are 
small. This is illustrated in Fig. |3| which shows the same 
string models shown in Fig. [3 with our Advanced LIGO 
amplitude estimate of ^50% = 10~^^ s~^/'^. 

To summarise, we find the chances of detecting "clas- 
sic" strings to be significantly smaller than previous es- 
timates suggest. Even Advanced LIGO only has a few 
percent chance of detecting "classic" strings at the mat- 
ter era peak (see the thin solid line around G/i ~ 10~^ in 
Fig. though it has a good chance of detecting cosmic 
superstrings and cosmic strings with small loops as show 
in Fig. 13 Initial LIGO requires the small reconnection 
probability of cosmic superstrings and/or small loops to 
attain a reasonable chance of detection. It should be 
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pointed out that the "classic" string model may well be 
incorrect. If loop sizes are given by gravitational back 
reaction then a reasonable guess for the spectrum of per- 
turbations on long strings leads to n = 3/2. This value, 
as shown on Fig. [21 has a substantially larger chance of 
detection by Initial LIGO. The size of loops, however, 
may not determined by gravitational back-reaction after 
all. 



B. Large loops 

Recent simulations 0, |^ suggest that loops are 
produced at sizes unrelated to the gravitational back- 
reaction scale. One of the simulation groups 19] finds a 
power law for the loop distribution, while the other two 
groups [23, 1^ find loops produced at a fixed fraction 
of the horizon, with loop production functions peaking 
around a « 10"'^ ^Q}, and the significantly larger a w 0.1 
l2U|. It should be pointed out that the first two results 
|l9ll2C| are expanding universe simulations, whereas the 
results of the third group |^ come from simulations in 
Minkowski space. 

Following formation, the length of loops shrinks due to 
gravitational wave emission according to , 

l{t) = k-TGii{t-U), (66) 

where 1^ — at^, is the initial length, and ti is the time of 
formation of the loop. The length goes to zero at time 

Loops are long-lived when tf ^ ti, i.e. when a/(rG/i) ^ 
1. For aw 0.1, using F = 50 the lifetime of loops is long 
provided Gfi ^ 2 x 10^'^, which covers the entire range 
of astrophysically interesting values of Gfi. On the other 
hand, if we take a « 5 x 10"'', then loops are long-lived 
only when Gfi < 10^^. 

If the size of loops is given by gravitational back re- 
action, then a is given by Eq. 1)5711 . and provided n > 1 
all loops are short-lived. This means we can use t) oc 
6(1 — at) (as we have so far) because the loop distribution 
is dominated by the loops that just formed. 

If loops are long-lived, the distribution can be calcu- 
lated if a scaling process is assumed (see §9.3.3 and 
§10.1.2). In the radiation era it is 

n{l,t) = Xr^"^^^(' + ^G'A^^)"^/^ 

I <at, t < tcq (68) 

where Xr ~ 0.4^0;^/^, and C is a parameter related to 
the correlation length of the network found in numerical 
simulations of radiation era evolution to be about 15 (see 
Table 10.1 in 0). The upper bound on the length arises 
because no loops are formed with sizes larger than at. 

In the matter era the distribution has two components, 
loops formed in the matter era and survivors from the ra- 
diation era. Loops formed in the matter era have lengths 
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FIG. 4: Plot of the rate of bursts 7 (top panel), and r\, 
the probability of having at least least one event in our data 
set with amplitude larger than ^50% in a year of observation 
(bottom panel), as a function of G/i for various loops distri- 
butions. For all curves we have set P = 50, /, = 75 Hz, 
c = p = 1, and the ignorance constants (?i = (72 = 1. As 
a reference we again show 7 and r\ computed using the loop 
distribution from Eq. 15611 . according to Eq. 159L in the A- 
universe with a — PG/i using the solid curve. The remaining 
curves have been computed through Eqs. 15511 and 13811 with 
^50% = 10"^" s"^''^. The dashed curves show 7 and rj com- 
puted using the loop distribution of Eqs. I68II . 169L and 1701 . 
The dashed-dot curves show 7 and rj computed with Eq. 1)71^ . 

distributed according to, 

ni{l,t) = Xmt-^{l + '^Gfity^ 

atcq — TG^{t — tcq) < I < at, t > ieq(69) 

with Xm « 0.12c, with C ~ 4 (see Table 10.1 in 12]). 
The lower bound on the length is due to the fact that 
the smallest loops present in the matter era started with 
a length ateq when they were formed and their lengths 
have since decreased due to gravitational wave emission. 
Additionally there are loops formed in the radiation era 
that survive into the matter era. Their lengths are dis- 
tributed according to, 

n2{l,t) = Xrtl^^t'Hl + rGfit)-^^^, 

Z < aicq - rG/i(t - <cq), t>tcq, (70) 

where the upper bound on the length comes from the fact 
that the largest loops formed in the radiation era had a 
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Eq. (71),|3=5/2, B=0.15 
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FIG. 5: Same as Fig. 0|but using our Advanced LIGO am- 
plitude estimate ^59% = 10^^""^ s^"'^''''. 

size atgq but have since shrunk due to gravitational wave 
emission. 

The simulations in [l9l | find a power law for the loop 
distribution, n{l) cx . If all loops produced are long- 
lived, meaning no loops are produced with sizes below 
TGfit, the distribution we expect is 

n{l,t) ^ Bt'^-^il + TGfit)-'^. (71) 

The most recent fits to their loop distribution find /3 « 
5/2, in both the radiation and matter eras, and B « 
0.1, 0.2 in the matter and radiation eras respectively [ssj . 

We can use the results of Sec. IV to compute the rate 
of bursts for these loop distributions. Figure 01 shows 
the burst rate 7, and the probability 77 of having at least 
least one event in our data set with amplitude larger 
than ^50% in a year of observation as a function of G/i 
for all the above distributions. For all curves we have 
set r = 50, /, = 75 Hz, c = p = 1, and the igno- 
rance constants gi =32 = 1- We have evaluated the 
cosmological functions in the A universe. As a refer- 
ence, we again show 7 and t] computed using the loop 
distribution from Eq. (|56l) . according to Eq. (|59f) . with 
a — TGfi using the solid curves. The remaining curves 
have been computed through Eqs. l(5B|) and with 
^50% = 10~^° s~^/^. The dashed curves show 7 and rj 
computed using the loop distribution of Eqs. Ht)8|) . (|69|l . 
and (|70|) with a = 0.1. The dashed-dot curves show 



7 and T] computed with Eq. H71|l . where we have taken 
/3 = 5/2, and B = 0.15 as an approximation for both the 
radiation and matter eras. Figure [S] shows the same loop 
distributions as Fig. 01 but using our Advanced LIGO 
amplitude estimate ^50% — 10^^^ s~^/^. 

The loop distributions shown here lead to a significant 
enhancement in the rate. Note that we have not included 
the enhancement in the rate of bursts for the case of 
cosmic superstrings p < 1. Given the range of results, 
it is important to determine whether the loop sizes at 
formation are determined by gravitational back-reaction, 
or whether they are large and we require a revised loop 
distribution. 



VI. RESULTS II: CONSTRAINTS 

If no events can be positively identified in a search, we 
can use the results of Sec. Ill to constrain the parameter 
space of theories that lead to the production of cosmic 
strings. Unfortunately, as we have mentioned, consider- 
able uncertainties remain in models of cosmic string evo- 
lution. Nevertheless, we can place constraints that are 
correct in the context of a particular string model. Here 
we will illustrate the procedure for the loop distribution, 
Eq. and the resulting rate Eq. 

We would like to absorb the ignorance constants gi 
and 52 into the parameters of the model. The two igno- 
rance constants enter the expression for the rate Eq. H59|) 
in three ways. First, they affect the upper limit of the 
integral Eq. Ht)l|) through Eq. Ht)0|l . which is, 

2/3/ \ A u-1/3 

{l + zy/^ifr{z) giGfiaV^' ^ > 

Secondly, they enter through 9m in the theta function 
cutoff of the rate, 

em = (.g2(l + z)f,aH^'Mz))~'/^ (73) 

and finally, the rate itself (Eq. (|59|l ) is proportional to 

-2/3 
.92 ■ 

If we write a = e(rG/i)", and substitute into Eqs. I|73|) 
and (EH), we can simultaneously absorb gi and 172 into 
new variables e{e, gi, g2,n) and X{Gijl, gi, g2)- Looking 
at Eq. (|73|) we can write an identity for the ignorance 
constants and the variables we want to absorb them into, 

g2e{TG^lr - e{TXr. (74) 

Similarly, looking at the denominator of the right hand 
side of Eq. (|72|l we write, 

<?iGAi[e(rG/x)'f = X[e(rA:)"]2/3 (75) 

These equations can be simultaneously solved to give, 

G^Ji = g^^gTx, (76) 
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Region ruled out 
at the 99.7% level 



Allowed Region 



FIG. 6: Contour of excluded regions at various confidences. 
In the cosmic string model we have set e = 1, and n = 3/2, 
and allowed X and /3 to vary. For the loudest event we have 
used our Initial LIGO estimate for the amplitude ^59% = 
10-2° s-i/^ 



and 



„ -2n/3-l- 



(77) 



If we replace G^i by Eq. (jTHJl and e by Eq. ^77^, in Eq. ifT^ 

we obtain, 

- ((1 + ^)^e(^X)"iy-V^(^))-'/^ (78) 
and if we make the same replacements into Eq. (|72|l . 



2/3/ N 



AH, 



-1/3 



(1 + z)i/3^^(^) X{e{TXy'Y/^'' 



(79) 



namely, we obtain functions of e and X only. 

Finally, if we replace by Eq. and e by Eq. (|77|l . 
in our expression for the rate, Eq. H59|l . we see that we 
can absorb the remaining factors by defining a quantity. 



« 1/3 C 

P = 52 51-- 



(80) 



This yields a rate. 



X 6(1 -0,„), 



(81) 



where 9m given by Eq. (|78|l . The parameters that we 
constrain are now X, n, e, and /3. 

Figure illustrates the application of this procedure. 
We show a contour plot constructed by computing the 
rate and comparing it to 768% ~ 1.14/r, 795% ~ 3.00/T 
and 799.7% ~ 5.81/T. In the cosmic string model we have 
set e = 1, and n = 3/2, and allowed X and j3 to vary. In 
place of the loudest event we have used our Initial LIGO 
estimate for the amplitude ^50% = lO"^" s"^/'^. If we 
set the ignorance constants to unity, and c = 1, we see 



that for in the range 4 x 10 ^ — 10 ^, values of the 
reconnection probability around p — are ruled out 
at the 99.7% level. 

It is interesting to note that for cosmic superstrings a 
bound on the reconnection probability might be turned 
into a bound on the string coupling and/or a bound on 
the size of extra dimensions. 



VII. SUMMARY AND CONCLUSIONS 

We begun this work by considering the data analysis 
infrastructure necessary to perform a search for grav- 
itational waves from cosmic (super)string cusps. The 
optimal method to use in such a search is matched- 
filtering. We have described the statistical properties of 
the matched-filter output, a method for template bank 
construction, and an efficient algorithm to compute the 
matched-filter. When data from multiple interferometers 
is available, we discussed consistency checks that can be 
used to greatly reduce the false alarm rate. 

The relevant output of the matched-filter is the esti- 
mated signal amplitude. We have shown that in a search, 
the event with the largest amplitude can be used to to set 
upper limits on the rate via the loudest event method. 
The upper limit depends on the cosmological rate of 
events and the efficiency, which can be computed using 
simulated signal injections into the data stream. The in- 
jections should account for certain properties of the popu- 
lation, such as the frequency distribution {dR/df oc f~^^^ 
for cusps), and source sky locations. Wc have briefly dis- 
cussed the related issue of sensitivity, which also depends 
on the search detection efficiency, as well as background 
estimation. We have made a single estimate for the am- 
plitude of detectable and loudest events using the Initial 
LIGO design curve and an SNR threshold pth — 4, which 
we believe is a reasonable operating point for a pipeline. 
Assuming that there are no significant inefficiencies in the 
pipeline, that the rate of false alarms is sufficiently low, 
and that no signals or instrumental glitches are present 
in the final trigger set, we expect this amplitude to be 
^50% ~ lO-^^s"^/^. Advanced LIGO is expected to be 
an order of magnitude more sensitive then Initial LIGO, 
so for Advanced LIGO we expect ^59% ~ lO-^^s"^/"^. 

We have computed the rate of bursts in the case of an 
arbitrary cosmic string loop distribution, using a generic 
cosmology. Motivated by the data analysis, we cast the 
rate as a function of the amplitude. We have applied 
this formalism in the case of a fiat universe with matter, 
radiation and a cosmological constant, to the loop distri- 
bution used in the estimates of 0,0, For the "clas- 
sic" (field theoretic strings with loops of size I = YGjit) 
model of cosmic strings considered in JJ„ J^J, we find 
substantially lower event rates than their estimates sug- 
gest. The bulk of the difference arises from our estimate 
of a detectable amplitude. The effect of the amplitude 
estimate on the rate is discussed in detail in Appendix 
B. The remaining discrepancy arises from differences in 
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the cosmology (the cosmological model in 0, 0, ^| in- 
cluded only matter and radiation), as well as factors of 
0(1) that were dropped in the previous estimates. The 
net effect is that the chances of seeing an event from 
"classic" strings using Initial LIGO data have dropped 
from order unity to about 10^'^ at the local maximum 
of the rate (which comes from bursts produced in the 
matter era). This result is illustrated in Fig. ^ 

The "classic" string model may well be incorrect so 
these results are not necessarily discouraging. Two new 
results indicate that this is the case. 

Recently it was realised that damping of perturbations 
propagating on strings due to gravitational wave emission 
is not as efficient as previously thought ■ As a conse- 
quence, the size of the small-scale structure is sensitive 
the spectrum of perturbations present on the strings, and 
can be much smaller than the canonical value TGfit |32j | . 
If the size of loops is given by gravitational back-reaction, 
a reasonable guess for the spectrum of perturbations on 
long strings leads to 

a ^ {TG^lf^ , (82) 

which in turn could lead to detectable events by Ad- 
vanced LIGO, or Initial LIGO if we are dealing with 
cosmic superstrings. This is illustrated in Figs. [21 and 

m 

Even more recently, simulations 0, |^ suggest 
that loops are produced at sizes unrelated to the gravi- 
tational back-reaction scale, a < 1. In this case loops of 
many different sizes are present at any given time (be- 
cause they are very long-lived) , and the use of a revised 
loop distribution becomes necessary. The particulars of 
the distribution are currently under debate, but all dis- 
tributions currently considered lead to an enhanced rate 
of bursts relative to the "classic" model. Results for the 
rate for various loop distributions are shown in Figs. 01 
and [SI 

Finally, we have have shown how the parameter space 
of theories that lead to the production of cosmic strings 
might be constrained in the absence of a detection. Re- 
sults for a model of cosmic superstrings, where the loop 
size is given by Eq. (|82|l . are shown in Fig. [Bl Initial 
LIGO may yield interesting constraints on the reconnec- 
tion probability for some range of string tensions. In- 
triguingly, a bound on the reconnection probability might 
be turned into a bound on the string coupling and/or a 
bound on the size of extra dimensions, in the case of cos- 
mic superstrings. 

Wish list 

Unfortunately, there remain considerable uncertainties 
in models of cosmic string evolution. While we can es- 
timate the detectability, and place constraints, that are 
correct in the context of a particular string model, the 
current parameter space allows for a wide spectrum of 
burst rates in the interesting range of string tensions. 



We would like to finish by posing a set of questions to 
the various cosmic string simulation groups which, from 
the perspective of gravitational wave detection, would 
greatly improve predictability. 

They are: 

• What is the size of cosmic string loops? Are they 
large when they are formed, so we need to consider 
a loop distribution? If so, what is that distribu- 
tion? Is their size instead given by gravitational 
back-reaction? If so, what is the spectrum of per- 
turbations on long strings? 

• What is the number of cusps per loop oscillation? 
Is it independent of the loop size? 

• What is the size of cusps? In particular, what frac- 
tion of the loop length I is involved in the cusp? 

• What are the effects of low reconnection probabil- 
ity? In particular, what is the enhancement in the 
loop density that results? 
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APPENDIX A: COSMOLOGICAL FUNCTIONS 

To derive exact expressions for ^pt{z), (fr{z), and 
ifv{z)^ we begin with the evolution of the Hubble func- 
tion. The Hubble function is given by 

H{z) - H^Kz), (Al) 
with, for a flat universe, 

h{z) = (a„(l + zf + a,(l + z)4 + nt,f'^ . (A2) 

Here, Vli ~ pi{z — Q)/pc{z = 0) is the present energy 
density of the i'th component, relative to the critical den- 
sity. The subscripts m, r, and A stand for matter (dark 
and baryonic), radiation, and the cosmological constant, 
respectively. Since we assume the universe flat, Eil^ = 1. 
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We will use the set of cosmological parameters in 
[3^. which provide a good fit to recent cosmological 
data. The precise values of the parameters are not crit- 
ical, but we include them here for clarity. They are, 
Hq = 73 km s-^Mpc-i = 2.4 x IQ-^^ s'^, = 0.25, 
VLr = 4.6 X 10"^, and JIa = 1 - fi™ - fi^. 

To compute the relation between the time and the red- 
shift, we use the fact that dz/dt = — (1 -t- z)H ^ and write 



t = dt' = 



dz' 



,0 J, {l + z')H{z'y 

The dimensionless function (pt{z) of Eq. H44|l is thus, 

dz' 



(A3) 



'^t{z) 



{l + z')h{z')' 



(A4) 



In order to compute the amplitude distance as a 
function of the rcdshift, we consider null geodesies in an 
FRW universe. We can take the polar and azimuthal 
coordinates to be constant and let the radial coordinate, 
r, vary. In this case we have that dr/dt = —(1 + z), and 
thus dr/dz = 1/H. So we write 



r = / dr' - , , , . 
/o Jo Hiz') 



dz' 



(A5) 



We can therefore express the dimensionless function of 
Eq. (071), fr{z), as 



ipriz) 



dz' 
H^' 



(A6) 



Finally, we would like to derive an expression for the 
differential volume as a function of the redshift, dV/dz. 
The differential volume element is given by 

dV = a^it^dr sin OdOdcj) , 

where a{t) is the scale factor. Integrating over the polar 
and azimuthal coordinates and using dr/dz gives 



dV = Ana-^ity^dr 



{l + zfH{z) 



dz . (A7) 



Using Eqs. I|A1(I and ljA5|l gives the dimensionless func- 
tion of Eq. 



ipv{z) = 



{i + zYh[zy 



(A8) 



APPENDIX B: APPROXIMATE ANALYTIC 
EXPRESSION FOR THE RATE 

In this appendix we compute an approximate expres- 
sion for the rate as a function of the amplitude can be 
obtained using the interpolating functions introduced in 
[Hill El], Eqs. ins, inSl), and ini- The result is use- 
ful to understand the qualitative behaviour of the rate 
curves shown in Figs. ^ 121 and|31 



Using Eqs. (|46|l and H63() we can write the amplitude 
of an burst from a loop of length I at a redshift z as 



A 



G/i^^/3(l + z)2/3 
tr)Z 



(Bl) 



We take the size of the feature that produces the cusp 
to be the typical size of loops, l{z) = aH^^ipt{z), with 
^Pt{z) given by Eq. H62I) . so that, 

A ^ GtJia^'h^^'\-\l + z)-i/3(l + z/z,,)-^'\ (B2) 
We define a dimensionless amplitude a, 
A 



Depending on the value of z, the function a(z) has three 
different asymptotic behaviours, 

{z-i z<l 
Z-4/3 K 2 < (B4) 
Zef 2-^/3 Z > Zeq 

The regime where z <C 1 corresponds to cusp events 
occurring nearby, the regime where 1 <C <C z^q cor- 
responds to matter-era events, and the regime where 
z ^ Zeq corresponds to radiation era events. 
At z = 1, we define 

ai = a{z = 1) « 2^1/3 - 1, 

and at z = Zeq, 

aeq = a{z = Zeq) « z-//='2-i/3 ^ z-g^/^ai « 4 x 10"^ 

This means we can write, 

(B5) 



^e<3 "eq 



The regime where z ^ 1, corresponds to a ^ 1; the 
regime where 1 <C z ^ Zg^, corresponds to a^q ^ a <C 1; 
and the regime where z 3> z^q, corresponds to a <C a^q. 
This means that we can write, 

z(a)^<( a-3/4 ae5«a«l (B6) 

a>l 

We can write this in terms of the interpolating function. 



/ X -3/20 

.(a) ~ a-3/2"a-3/^ [l + j (1 + a) 



-1/4 



(B7) 



Using Eqs. Ht)2|l and 1)64(1 . we can write the rate of 
events as a function of the redshift, Eq. (|59|) . as, 



dR 

dz 



6z2(l^-z)-^/6(l^-z/ze,)"/^ (B8) 
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where b is defined as 



(B9) 



and the theta-function cutoff has been ignored for con- 
venience. In the three regimes considered above the rate 
is given by, 

{z^dz z < 1 

z^/^dz K z < (BIO) 
Z-,"/'z8/3dz Z » Ze, 

In the regime where z <JC 1, corresponds to a 3> 1, and 
z ^ a^^. This means, 



Since, 



we have that 



dz 

dR ~ bz'^dz = ba^^ —da. 

da 



dz 
da 



dR ^ —ba '^da, for a 3> 1. 



(Bll) 



(B12) 



(B13) 



The regime where 1 ^ z ^ z^q, corresponds to aeq <C 
a 1, and z ^ a^'^/*. This means, 

d7 

dR ~ fez^/^dz = ba-'^^^—da. (B14) 
da 



Since, 



dz _3 y/4 



(B15) 



we have that 

dR ^ba'^^/^da, for < a < 1. (B16) 

Finally, the regime where z ^ z^q, corresponds to a <C 
tteq, and z ~ (aeq/ai)~^/2°a~^/^. Thus, using Eq. I|B5|1 . 

di? ~ 6z-"/'^'/'rf^ = &a39/40a-«/5^da. (B17) 
^ ^ da 



Since, 



dz 3 



3/20^-8/5 



(B15 



da 5 
we have that 

di? ~ -f 6a^i^/^Vi6/5da, for a<^a,q. (B19) 

Summarising, 



di?(a) 
da 



-b X 



33/40 -16/5 



-19/8 ae, « a « 1 (B20) 
i-^ a > 1 



Eq. (|B20|) can be integrated, to give the rate of events 
with reduced amplitude greater than a. 



5x 



33/40 _ 
aeq a 



11/5 



a < aeq 
aeq <C a <C 1 
a > 1 



(B21) 



We can determine the functional dependence of the 
rate for the simple case when a — TGfj,, which we can 
then compare with the dashed and dashed-dot curves of 
Fig.n The factor of b, as defined by Eq. (|B9p . contains a 
factor of {Gfj,)~^ as well as a factor of {G^)~^^^ through 
its dependence on a. So we take b oc [Gijl)-^''^. The 
dimensionless amplitude a, as defined by Eq. (jBSp also 
contains a factor of {G^)~^ as well as a factor of (G/i)"^/'^ 
through its dependence on a, so that a cx (G/i)"^/'^. 

In Eq. (jB21|) , (as we have mentioned) the regime where 
a <C aeq maps into the radiation era, and in this case the 
rate, 



R cx G/i. 



(B22) 



The regime where aeq <C a <C 1 maps into the matter era 
when the redshift z 1, and the rate in this case is, 



Ro,{G^lr^^'^\ 



(B23) 



The final regime when a ^ 1 corresponds to bursts that 
are coming from close by (z ^ 1), and the rate. 



i?oc (G^)^/^ 



(B24) 



So we immediately see that the first steep rise in the rate 
as a function of G/i (from left to right) of the dashed and 
dashed-dot curves of Fig. Q corresponds to bursts that 
are coming from small redshifts, i.e. from z <^ 1. The 
slight decrease in the rate comes from bursts produced at 
large redshifts, but still in the matter era, i.e. 1 ^ z ^ 
Zgq, and the final increase in the rate comes from bursts 
produced in the radiation era, z ^ Zeq. 

Eq. ljB21|l also makes it easy to understand the lower 



burst event rates we find relative to the previous esti- 
mates of Damour and Vilenkin ll m m (see the the 
dashed and dashed-dot curves of Fig. . They compared 
the strain produced by cosmic string burst events at a 
rate of 1 per year to a noise induced SNR 1 event, given 
the Initial LIGO design noise curve (see, for example, the 
dashed lines in Fig. 1 of ^3)- I'^ ou'" treatment, the am- 
plitude that corresponds to is A w lO^^^s^^/'^. To make 
our amplitude estimate we have chosen an SNR threshold 
of 4, which, as proposed in is a reasonable operating 
point for a data analysis pipeline, and have included the 
effects of the antenna pattern of the instrument, which 
averages to a factor of VS. The combined effect is to 
increase the amplitude estimate by about a factor of 10, 
to make it A « lO-^^s-^/^. 



Looking at Eq. (|B21|) . we see that an increase in the 
amplitude of a factor of 10, results in a decrease of a 
factor of lO'^ in the rate of nearby bursts (z ^ 1), a 
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decrease by a factor of about 24 in the rate of bursts by about a factor of 160 in the rate of bursts produced 
produced in the matter era {1 <^ z -^eg); ^md a decrease in the radiation era (z ^ z^g). 
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